Direct measurements of anisotropic energy transfers in a rotating turbulence 

experiment 
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We investigate experimentally the influence of a background rotation on the energy transfers in 
decaying grid turbulence. The anisotropic energy flux density, F(r) = (^u(^u)^), where is the 
vector velocity increment over separation r, is determined for the first time using Particle Image 
Velocimetry. We show that rotation induces an anisotropy of the energy flux V • F, which leads to 
an anisotropy growth of the energy distribution ^(r) = ((^u)^), in agreement with the Karman- 
Howarth-Monin equation. Surprisingly, our results prove that this anisotropy growth is essentially 
driven by a nearly radial, but orientation-dependent, energy flux density F(r). 



The energy cascade from large to small scales, and the 
associated Kolmogorov 4/5th law, are recognized as the 
most fundamental results of homogeneous and isotropic 
turbulence [HE]. In the presence of a background rota- 
tion, a situation which is relevant for most geophysical 
and astrophysical flows, the scale-to-scale energy trans- 
fers are modified by the Coriolis force, yielding a grad- 
ual columnar structuring of turbulence along the rota- 
tion axis The Taylor-Proudman theorem is of- 
ten invoked, however improperly, to justify the resulting 
quasi-2D nature of turbulence under rotation. Indeed, 
this theorem is a purely linear result, which applies only 
in the limit of zero Rossby number (i.e. infinite rotation 
rate), and is therefore incompatible with turbulence; it 
cannot describe the anisotropic energy transfers responsi- 
ble for the non-trivial organization of rotating turbulence 
which are a subtle non-linear effect taking place only at 
non-zero Rossby number. To date, no direct evidence for 
these anisotropic energy transfers towards the 2D state 
in the physical space has been obtained. In this Let- 
ter, we report for the first time direct measurements of 
the physical-space energy transfers in decaying rotating 
turbulence using Particle Image Velocimetry (PIV), and 
provide new insight into the anisotropy growth of turbu- 
lence at finite, and hence geophysically relevant, Rossby 
number. 

If homogeneity (but not necessarily isotropy) holds, the 
energy distribution and energy flux density in the space 
of separations r are described by the fields 



and F(r,t) = ((^u((^u)^), (1) 



where u(x, t) is the turbulent velocity, = u(x + r, — 
u(x, t) is the velocity vector increment over r (Fig. [T]), 
and (•) denotes spatial and ensemble averages. These key 
quantities satisfy the Karman-Howarth-Monin (KHM) 
equation (TJ |8] , which describes the evolution of the en- 
ergy distribution in the space of separations. 
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where i?(r,t) = (u(x, t) • u(x + r, t)) = (u^) - E(r, t)/2 is 
the two-point velocity correlation and v the kinematic 




FIG. 1: (a) Experimental setup. The water- filled tank is ro- 
tating at < n < 1.68 rad s~^. The grid is towed from the 
bottom to the top, and PIV measurements are performed in 
the vertical plane (x, z) in the rotating frame during the tur- 
bulence decay, (b) Definition of the vector velocity increment 
= u(x + r) — u(x). 

viscosity. Importantly, this equation is still valid for 
homogeneous anisotropic turbulence [9 , and in partic- 
ular for axisymmetric turbulence in a rotating frame 
(here axisymmetry is to be understood in the statisti- 
cal sense, with respect to r). For stationary (forced) 
turbulence, this equation reduces to V • F = — 4e in 
the inert ial range, where e stands for the rates of in- 
jected and dissipated energy. In the isotropic case, this 
constant-flux relation yields a purely radial flux density, 
F(r) = — (4/3)er, describing the usual energy cascade 
from large to small scales. This result is actually identi- 
cal to the celebrated Kolmogorov's 4/5th law, classically 
expressed in terms of the 3rd order longitudinal structure 
function, {5u\) = — (4/5)er, where Sul = Su-r/r is the 
longitudinal velocity increment. 

In decaying rotating homogeneous turbulence, Eq. ([2| 
shows that, starting from an isotropic initial energy dis- 
tribution £^(r, 0), an anisotropy growth in E(Y^t) is ex- 
pected if an anisotropic energy flux V-F is induced by the 
Coriolis force. However, the flux density F(r) itself has 
never been measured, and its precise form, which reveals 
the fundamental action of rotation on turbulence, is so 
far unknown. The only experimental attempts to charac- 
terize the energy transfers in rotating turbulence were re- 
stricted to measurements of {5u\) in the plane normal to 
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the rotation axis [TOl [TT] , hence ignoring the anisotropic 
nature of those transfers. Recent theoretical efforts have 
been made to generahze the 4/5th law, assuming weak 
anisotropy [12 , or considering the full anisotropic prob- 
lem but restricted to the stationary case 

Experiments.- The experimental setup is similar to 
the one described in Ref. [13] , and is briefly recalled here 
(Fig. [l^). Turbulence is generated by towing a square 
grid at a velocity Vg = 1.0 m from the bottom to the 
top of a cubic glass tank, of side 60 cm, filled with 52 cm 
of water. The grid consists in 8 mm thick bars with a 
mesh size M = 40 mm. The whole setup is mounted on 
a precision rotating turntable of 2 m in diameter. Runs 
for three rotation rates, Q = 0.42,0.84 and 1.68 rad s~^ 
(4, 8 and 16 rpm), as well as a reference run without ro- 
tation, have been carried out. The initial Reynolds num- 
ber based on the grid mesh is Rcg = VgM/v = 40 000, 
and the initial Rossby number Rog = Vg/2QM ranges 
from 7.4 to 30, indicating that the flow in the close wake 
of the grid is fully turbulent and weakly affected by ro- 
tation. During the turbulence decay, the instantaneous 
Rossby number, Ro{t) = (u^)^/^/2r^M, decreases with 
time down to 10~^, spanning a range in which influence 
of rotation is expected. An important concern about grid 
turbulence experiments in a confined rotating domain is 
the excitation of reproducible inertial modes [14^. Here, 
we use the modified grid introduced in Ref. [13], which 
was shown to significantly reduce the generation of these 
modes. Consequently, turbulence can be considered here 
as almost freely decaying and homogeneous, a necessary 
condition for the validity of the KHM equation ([2|. 

Velocity measurements are performed in the rotating 
frame using a corotating PIV system. Two velocity com- 
ponents {ux^Uz) are measured, in a vertical 16 x 16 cm^ 
field of view, where z is the rotation axis. During the 
decay of turbulence, 60 image pairs are acquired by a 
double-frame 2048^ pixels camera, at a rate of 1 pair per 
second. The PIV resolution, 1.3 mm, is sufficient to re- 
solve the inertial range but fails to resolve the dissipative 
scale (the Kolmogorov scale is of the order of 0.2 mm 
right after the grid translation [11 ). 

Only surrog ates of £;(r) and F(r) ^ can be computed 
from the measured 2D velocity fields. These surrogate 
quantities are defined as 

^(r) = {5ul + 5ul):,,,, F(r) = {5n{5ul + 5ul))^,,, (3) 

where the spatial average is restricted to the measure- 
ment plane, and r = Vx^x^'f^z^z- For each time after the 
grid translation, these quantities are computed for all 
separations r in the PIV field of view, and are ensemble- 
averaged over 600 realizations of the turbulence decay. 
The fields ^(r) and F(r) are remapped on a spherical 
coordinate system (r, ^,0), where r = |r|, and is the 
polar angle between Bz and r; the invariance with respect 
to the (non-measured) azimuthal angle (j) is assumed by 
axisymmetry. Although relations between the surrogates 
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FIG. 2: Energy distribution ^(r) at time tVg/M = 400 after 
the grid translation, for (a) Q = 0, and (b) Q = 1.68 rad s~^ 
(16 rpm). (c) Horizontal-to- vertical energy ratio as a function 
of time at scale r = 10 mm for various Q; o: additional curve 
at r = 30 mm for Q = 16 rpm. Stars indicate integer numbers 
of tank rotations. 

([3| and the exact 3-components quantities ([T]) can be de- 
rived for isotropic turbulence, no general relation holds 
in the anisotropic case^so we do not apply any correc- 
tion weight in E and F. Since only the surrogates are 
considered in this paper, we simply drop the tildes ^ in 
the following. 

The convergence of the statistics from experimental 
measurements is very delicate to achieve, in particular 
for the computation of F(r), which is a 3rd order mo- 
ment of a zero- mean velocity increment. We found that, 
using a set of 600 realizations of the turbulence decay, a 
convergence better than 5% at small scales, and of the 
order of 20% at scales r M, could be achieved for F(r). 
The convergence for E{r) is better than 1% for all scales 
up to r M. 

Energy distribution.- The map of energy distribution 
E{r) for separations r in the vertical plane is plotted in 
Fig. [2j at a time tVg/M = 400 after the grid translation. 
The iso-£^ curves are found nearly circular for = (Fig. 
|2|i), showing the good level of isotropy of our grid tur- 
bulence without rotation. On the other hand, they are 
highly anisotropic at the same time for 1] = 16 rpm (cor- 
responding to 4.3 tank rotations), with a strong depletion 
of E{r) along the rotation axis z (Fig. ^p)- The deple- 
tion of E{t) corresponds to an enhanced velocity correla- 
tion R{t) along the rotation axis, reflecting the classical 
trend towards a 2D flow invariant along z. Importantly, 
an isotropic energy distribution is found in the 3 rotat- 
ing cases just after the grid translation, as demonstrated 
in Fig. |2|c) where the time evolution of the horizontal- 
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FIG. 3: Energy flux density F(r) in the non-rotating case, 
at time tVg/M = 400 after tfie grid translation, (a) Flux 
density magnitude |F|. (b) Raw vector field F. (c) Energy 
flux V • F. (d) Deviation angle a(r), defined as sina(r) = 
ey • (cr X F)/|F|; iso-angle lines are separated by 5°. The 
dashed line in (c) shows the "crest line", following the local 
maximum of — V • F. 

to- vertical energy ratio E{6 — {))/E{0 = 7r/2) is plotted 
for an inertial-range separation r = 10 mm. This con- 
firms that the initial grid turbulence is isotropic even 
when Q ^ 0^ and that the subsequent anisotropy growth 
is a pure effect of the background rotation. Fig. [2|c) 
also shows that the anisotropy growth rate is essentially 
proportional to ft [Sl [I] • Interestingly, the anisotropy is 
found more pronounced at small scales, as shown by the 
lower anisotropy ratio plotted for r = 30 mm. It is worth 
noting that this stronger anisotropy at small scales is in 
contradiction with the naive assumption that large scales, 
having a slower dynamics, are more affected by rotation 
than the faster and supposedly still 3D small scales. 

Energy transfers: isotropic case.- We now turn to the 
energy flux density, and we first present in Fig. [sj^b) 
measurements of F(r) for Q = 0, at the same time 
tVg/M = 400. This vector field is found remarkably ra- 
dial, pointing towards the origin, giving direct evidence of 
the isotropic energy cascade in the physical space, from 
the large to the small scales, in the non-rotating case. 
Finer assessment of the isotropy of F can be achieved 
by introducing the following three scalar quantities: the 
deviation angle a{r) from the radial direction (Fig. [sji), 
the magnitude |F| (Fig. [3^), and the energy flux V • F 
(Fig. [sj^). The very weak angle measured for r < M, 
a{r) 2° ± 2°, confirms the almost purely radial nature 
of F. The isotropy of the flux density magnitude is not 
as good: the iso-|F| are nearly circular up to r 30 mm, 
but shows slight departure from isotropy at larger r, sug- 
gesting that this quantity is very sensitive to a residual 
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FIG. 4: Energy flux density F(r) in the rotating case (Q = 
16 rpm), at time tVg/M = 400 after the grid translation. 
Same layout as for Fig. [3] 

anisotropy of the large-scale flow. However, the iso-V • F 
remain remarkably circular up to r M, showing that 
the residual large-scale anisotropy has indeed a weak in- 
fluence on the energy flux for r < M. The energy flux 
V • F shows a broad negative minimum in an annular 
region spanning over r 5 — 20 mm, providing an indi- 
cation of the extent of the inertial range (we recall that, 
in the inertial range, V • F = — 4e), and decreases to zero 
both at small and large scales. 

Energy transfers: rotating case.- We consider now 
the energy transfers in the rotating case, shown in Fig. [4] 
at the same time tVg/M = 400. Interestingly, the flux 
density F is found to remain nearly radial for all sepa- 
rations, in qualitative agreement with recent predictions 
[9 , except at the smallest scales, for r < 10 mm, where a 
marked deflection towards the rotation axis is observed. 
Such horizontally tilted F is indeed consistent with an 
asymptotic 2D flow, for which F must be a strictly hori- 
zontal vector, function of the horizontal component of the 
separation only. This small-scale anisotropy is best ap- 
preciated from the map of the deviation angle a (Fig.jiji), 
showing a region of nonzero a at small scale only. Note 
that horizontally tilted F exists only for intermediate an- 
gle since axisymmetry requires a radial F for ^ = and 
7r/2. The 2D trend is remarkably weak in terms of the 
orientation of F(r) in the inertial range, compared to 
the strong anisotropy observed for the energy distribu- 
tion E{r) at comparable scales: a is only in the range 
— 10° in the inertial range, and increases up to 25° ± 5° 
for r ^ 0, with no significant dependence with Q. 

If we focus on the flux density magnitude |F|, which 
is essentially given by the radial component —Er = — F • 
e^, a clear anisotropy is now found at all scales. This 
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FIG. 5: (a) Spatio-temporal diagram of the horizontal-to- 
vertical energy flux ratio, V-F(^ = 7r/2)/V-F(^ = 0), showing 
the anisotropy growing from small to large scales, (b) Time 
evolution of the energy flux ratio at scale r = 10 mm for 
various Q. 

suggests that the anisotropy of the energy transfers is 
mostly driven by the ^-dependence of F^, and not by the 
growth of a nonzero polar component Fq = ¥ ■ eg. The 
maximum of |F| is systematically encountered near the 
rotation axis, at rather large scales, centered around SO- 
SO mm (outside the range shown in Fig. [4^). The local 
maximum of |F| on the horizontal axis is encountered at 
smaller scales, as evidenced by the crest line in Fig.|4ja). 

The flux map V • F (Fig. ^) shows an overall 
anisotropic structure similar to that of |F|, but essentially 
shifted towards smaller scales. The inert ial range, where 
the flux V • F is negative and approximately constant, 
becomes vertically elongated as time proceeds. Actually, 
although |F| is maximum along the rotation axis, it is 
spread over a wider range of scales, leading to a weaker 
flux V • F along z than along x, and hence a less intense 
vertical energy cascade. Here again, this is consistent 
with a 2D trend, which should yield a vanishing energy 
flux along the rotation axis. The horizontal-to- vertical 
flux ratio in Fig. |5] illustrates this vanishing vertical en- 
ergy cascade as time proceeds, an effect which is clearly 
enhanced as the rotation rate is increased. 

It must be noted that the spatial structure of the flux 
V • F is in good qualitative agreement with the KHM 
equation ([2|. Indeed, neglecting the viscous term, the 
vertically elongated region where V • F < induces a 
stronger reduction of the velocity correlation R along 
X than along resulting in a relative growth of the 
vertical correlation along and hence a vertical deple- 
tion of the energy distribution E = 2((u^) — R). We 
can conclude that the measured flux density F contains, 
through its divergence, a spatial structure consistent with 
the anisotropy growth of E observed in Fig|2j Interest- 



ingly, in line with the stronger anisotropy of E{t) found 
at smaller scales, the flux is also found more anisotropic 
at smaller scales. This is clearly demonstrated by the 
spatio-temporal diagram in Fig. [sja), showing that the 
anisotropy first appears at small scales, and then propa- 
gates towards larger scales as time proceeds. 

Conclusion.- We report the first direct measurements 
of the energy flux density F in the physical space in a 
decaying rotating turbulence experiment. Although the 
alternative description of the energy transfers in the spec- 
tral space is more natural for theory or numerics [2]-[4l [6] , 
the direct use of the KHM equation ([2| in the physical 
space, which is better suited for experiments, reveals here 
new and unexpected behaviors. The spatial structure of 
the measured energy distribution and energy flux V-F are 
found in good qualitative agreement with the KHM equa- 
tion which, to our knowledge, has never been assessed 
experimentally. Surprisingly, the anisotropy growth of 
the energy distribution is primarily driven by an almost 
radial, but orientation-dependent, flux density F, except 
at small scales where F shows a horizontal tilt, compat- 
ible with a trend towards a 2D state. It is also demon- 
strated that the anisotropy is paradoxically stronger at 
small scales, and propagates towards larger scales as time 
proceeds, an unexpected result which should motivate 
new theoretical efforts. 
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[1] 

[2] 

[3] 

[4] 
[5] 

[6] 

[7] 
[8] 

[9] 
[10] 

[11] 
[12] 

[13] 
[14] 



U. Frisch, Turbulence - The Legacy of A. N. Kolmogorov 
(Cambridge University Press, Cambridge, 1995). 
P. Sagaut and C. Cambon, Homogeneous Turbulence Dy- 
namics (Cambridge University Press, Cambridge, 2008). 
C. Cambon and L. Jacquin, J. Fluid Mech. 202, 295 
(1989). 

F. Waleffe, Phys. Fluids A 5, 677 (1993). 

P.J. Staplehurst, P.A. Davidson, and S.B. Dalziel, J. 

Fluid Mech. 598, 81 (2008). 

PD. Mininni and A. Pouquet, Phys. Fluids 22, 035106 
(2010). 

F. Moisy et al, J. Fluid Mech. 666, 5 (2011). 

A. S. Monin and A. M. Yaglom, Statistical Fluid Mechan- 
ics, vol. 2 (MIT Press, Cambridge, 1975). 
S. Galtier, Phys. Rev. E 80, 046301 (2009). 
C. N. Baroud et al, Phys. Rev. Lett. 88, 114501 (2002). 
C. Morize, F. Moisy, and M. Rabaud, Phys. Fluids 17, 
095105 (2005). 

S. Chakraborty and J. K. Bhattacharjee, Phys. Rev. E 
76, 036304 (2007). 

C. Lamriben et al, Phys. Fluids 23, 015102 (2011). 

G. P. Bewley et a/., Phys. Fluids 19, 071701 (2007). 



